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FREQUENCY  SHIFTS  INDUCED  IN  LASER  PULSES 
BY  PLASMA  DENSITY  VARIATIONS 

The  propagation  of  electromagnetic  radiation  through  plasmas  is  a  problem  of  gen¬ 
eral  interest  with  a  wide  variety  of  applications  ranging  from  communications^  to  laser 
driven  particle  accelerators.^’®  For  example,  the  results  of  plasma  simulation  studies  have 
recently  been  reported^’®  which  suggest  two  possible  methods  by  which  the  frequency 
of  an  electromagnetic  wave  may  be  upshifted.  In  the  first  method,  the  plasma  density 
through  which  the  EM  wave  is  propagating  is  suddenly  increased  in  time,^  while  the  sec¬ 
ond  method  utihzes  the  interaction  of  a  plasma  wave  (having  a  phase  velocity  near  the 
speed  of  light with  a  short  EM  pulse.®  Phenomena  such  as  these,  which  result  from 
variations  in  the  plasma  density,  may  offer  a  way  of  tuning  the  radiation  from  a  laser  or, 
alternatively,  they  may  describe  the  distortion  of  radio  signals  in  the  ionosphere. 

In  the  following,  a  general  analytic  theory  is  developed  from  first  principles  which 
describes  how  a  EM  radiation  pulse  is  affected  by  variations  in  the  plasma  density.  Specif¬ 
ically,  the  ID  wave  equation  is  used  to  derive  expressions  for  the  shifts  induced  in  the 
frequency  and  ia  the  wavenumber  by  arbitrary  plasma  density  variations  (both  in  space 
and  in  time).  In  particular,  for  variations  in  the  plasma  density  as  a  lunction  of  time,  it  is 
shown  that  the  wavenumber  of  the  EM  wave  remains  constant  while  shifts  are  induced  in 
the  frequency.  (This  is  to  be  contrasted  to  the  case  of  spatial  plasma  density  variations, 
for  which  the  frequency  remains  constant  and  the  wavenumber  shifts.)  For  the  case  of  a 
short  radiation  pulse  interacting  with  a  plasma  wave  with  finite  phase  velocity,  it  is  shown 
that  shifts  are  induced  in  both  the  wavenumber  and  in  the  frequency.  Maximum  frequency 
shifts  may  be  obtained  when  the  phase  velocity  of  the  plasma  wave  is  equal  to  the  group 
velocity  of  the  radiation  pulse  evaluated  at  the  ambient  plasma  density.  Furthermore,  it 
is  shown  that  for  small  frequency  shifts,  the  amplitude  of  the  vector  potential  mEuntains 
its  initial  profile  while  propagating  forward  at  the  ambient  group  velocity. 

The  ID  wave  equation  for  the  normalized  vector  potential  a.  =  eA/{mc^)  of  the  radi¬ 
ation  field  is  given  by 

(1) 

where  =  kpgn(z,t)/(y(z,t)no)-  Here  7(2, t)  is  the  relativistic  factor  associated  with  the 
motion  of  the  plasma  electrons,  n(z,i)  is  the  plasma  electron  density  and  kp^  =  Wpo/c®, 
where  Wpo  is  the  electron  plasma  frequency  in  the  ambient  density  tiq.  In  deriving  the 
above  equation,  conservation  of  canonical  momentum  was  used,  p±  =  cAx/c,  which  gives 
a  leans »erse  plasma  current  (in  the  fluid  limit)  of  J±  =  —enp±/(mj).  Throughout  the 
Manuscript  approved  September  8,  1989. 
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following  it  will  be  assumed  that  a?  <<  \  such  that  kp{z,t)  is  independent  of  a(z,t),  i.e., 
the  effects  of  the  radiation  field  on  the  plasma  wave  will  be  ignored. 

To  solve  the  above  wave  equation,  it  is  helpful  to  write  a(z,t)  =  b(z,t)exp{ikoz  — 
tu;ot),  where  b{z,t)  is  the  radiation  envelope  and  where  wo  and  ko  are  the  frequency  and 
wavenumber  of  the  radiation  in  the  ambient  plasma  (in  the  absence  of  a  plasma  wave) 
which  satisfy  the  dispersion  relation  =  c^kl  +  c^fcpQ.  Furthermore,  it  is  convenient  to 
introduce  a  change  of  variables  (  =  z  —  Vgt  and  t  =  t,  where  Vg  =  c^ko/u}o  is  the  group 
velocity  of  the  radiation  in  the  ambient  plasma.  The  wave  equation  is  then  given  by 


+  2^-^  4-  — — 
\  c^  dr  d(dT  7^ 


1  92  \ 
c2  dr^  ) 


b  ^  6kl{C,r)b{(,T), 


(2) 


v.-here  Sk^  =  k^  -  and  1/7=  =  1  -  v^'/c^. 


It  will  now  be  assumed  that  the  radiation  envelope  is  slowly  varying  compared 

to  the  radiation  frequency  u»o,  that  is  \db/dT\  <<  |tJo6|  and  \db/dC,\  <<  |u>o6/c|.  Assuming 
this,  the  second  order  derivatives  in  Eq.  (2)  may  be  neglected.  Equation  (2)  may  then  be 
solved  giving 


b{(,r)  =  6o«)exp 


tc* 


i: 


dr'SkliC^r') 


(3) 


[  2u^o 

The  condition  \db/dT\  <<  |u;ol>|  implies  \c^8kpl{2u)l)\  <<  1.  For  a  plasma  density  pertur 
bation  ^n,  this  condition  implies  |u;po5n/(2u;^no)|  <<  1.  The  condition  |96/5(|  <<  |u>o6/c 


implies  \  Jq  dT\d8k'^l dC,)c^ l{2uil)\  «  1  as  well  as  |(d6o/<fC)c/(i>oWo)|  <<  1.  The  first  of 
these  inequalities  generally  implies  that  the  wavenumber  shift  must  be  small  compared  to 
a;o/c,  while  the  second  inequality  indicates  that  the  initial  envelope  60(C)  must  be  slowly 
varying  compared  to  u^o  / c. 


The  total  phase  $(C,t)  of  the  radiation  field  may  be  identified  by  writing  a(C,T)  = 
|a(C,T)|  expi$(C,T).  It  is  then  possible  to  examine  the  evolution  of  the  frequency  as 
well  as  the  wavenumber  of  the  radiation  through  the  definitions  a;(C,7')  =  —d^jdi  = 
—{d/dr  —  Vgd/dC^)^  and  k{t^,T)  =  d^fdz  =  d^/d^.  Using  the  above  solution  for  6(C,t) 
gives 

=  dr’±6kl[(y),  (4) 

=  r<(T'^i*“(C,T').  (5) 

The  above  equations  are  veilid  for  arbitrary  variations  6kp((,T).  (Recall,  k^  =  fcpou/(uo7) 
and  the  effects  of  7  become  important  for  nonlinear  relativistic  plasma  waves. Provided 
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remains  slowly  varying  compared  to  a;o,  the  amplitude  of  the  normalized  vector 
potential  (a(  does  not  change.  That  is,  Eq.  (3)  indicates  that  |a)  =  |l>o(C)l  and,  hence,  the 
initial  envelope  of  the  vector  potential  is  simply  convected  forward  at  the  group  velocity 
Vs- 

To  illustrate  the  above  theory,  consider  a  plasma  density  variation  which  is  a  function 
only  of  space.  For  example,  consider  a  radiation  pulse  entering  a  plasma  <<  1) 

from  vacuum  with  a  plasma  density  profile  Sk^  =  8kp{z)  for  z  >  0  and  equal  to  zero 
for  2  <  0.  Assume  that  at  <  =  0,  the  radiation  pulse  extends  from  —L  <  z  <  0,  where 
woL/c  »  1.  Equations  (3)-(5)  indicate,  that  as  the  pulse  propagates,  the  frequency  and 
wavenumber  evolve  according  to  a>(2,<)  =  u>o  and  k{z,t)  =  ko  —  c^8k^{z)l{2vgU)o).  This 
is  in  agreement  with  the  well-known  result^  that  as  radiation  propagates  into  a  plasma 
with  spatial  density  variations,  the  frequency  remains  constant  whereas  the  wavenumber 
changes  such  that  the  dispersion  relation  =  c'^[k^  +  ^p)  remains  satisfied. 

On  the  other  hand,  consider  a  plasma  density  variation  which  is  a  function  only  of  time. 
For  example,  consider  a  radiation  pulse  (of  length  L  >>  c/u>o)  propagating  through  a  long, 
uniform  plasma  column  (where  <<  1)  in  which  the  density  is  temporally  changing, 

6k^  ~  Sk^{t).  Equations  (3)-(5)  indicate,  that  as  the  pulse  propagates,  the  frequency  and 
wavenumber  evolve  according  to  u>(2,t)  =  u>o  4-  c^6kp{z)/{2ijoo)  and  fc(2,t)  =  fco-  This  is 
in  agreement  with  the  simulations  of  Ref.  4  which  indicate  that  as  radiation  propagates 
through  a  plzisma  with  temporal  density  variations,  the  wavenumber  remains  constant 
whereas  the  frequency  changes  such  that  the  dispersion  relation  =  c^{k^  +  fcp)  remains 
satisfied. 


Consider  a  plasma  density  variation  which  is  a  function  both  of  time  and  of  space. 
For  example,  consider  the  case  in  which  a  plasma  wave  (with  phase  velocity  near 
is  used  to  upshift  the  frequency  of  a  laser  pulse,  as  suggested  by  the  simulations  of  Ref.  5. 
Assuming  a  plasma  wave  with  a  phase  velocity  Vp  such  that  the  plasma  wave  quantities 
are  functions  of  only  z  —  Vpt  implies  6kp{<^,T)  =  Skp{(^  —  Avr),  where  SkpiC  ~  Avr)  has 
the  form  of  a  periodic  oscillation  and  Av  z=Vp~Vg.  Defining  the  shift  in  frequency  and  in 
wavenumber  as  Aw  =  w(^,t)  —  w((,0)  and  Afc  =  k{(l,T)  —  fc(C,0)  gives 


Aw  = 


c"Vp 
2wo  Ar 


(6) 


Afc  = 


2wo  Av 


(T) 


The  above  equations  indicate  that  the  frequency  shift  will  be  maximum  for  the  case 


oist 


A'* 


Codes 

.'.v.M  p.nJ/or 
S'r’ec  lal 


□  □ 


Vp  =  Vg.  For  this  case  Aw  =  —{vgTc^l2u}^)dSkpld(,  and  Afc  =  — (c^r/2wo)d6fcp/dC-  Hence, 
Aw  is  a  linear  function  of  the  distance  that  the  pulse  has  propagated  through  the 
plasma.  In  such  a  way  jAw]  will  increase  until  it  becomes  sufficiently  large  so  that  the 
approximation  that  b((,T)  is  slowly  varying  is  no  longer  veilid.  For  6(C,t)  to  be  slowly 
varying  requires  |cAw/(vjWo)P  <<  1.  Consider  a  plasma  wave  with  Vp  =  Vg  such  that 
Skp(()  =  kpoSn(()/no,  where  in((^)  is  the  plasma  wave  density  perturbation  For  6n(C)  = 

Sno  sinfcpoC?  f 

Aw  <*>poSnoVgT 

-  = 1 - r^coskpoC-  (8) 

Wo  Wq  no  Apo 

Hence,  a  positive  frequency  shift  Aw  >  0  requires  the  laser  pulse  to  be  positioned  prop¬ 
erly  in  the  phase  of  the  plasma  wave  such  that  dSn/dC,  <  0.  This  is  illustrated  in  Fig. 
1.  A  convenient  method  for  producing  such  a  plasma  wave  may  be  the  laser  Wakefield 
accelerator,®  in  which  a  driving  laser  pulse  of  frequency  wq  is  used  to  generate  a  plasma 
wave  with  Vp  =  Vg. 

For  the  case  Vp  ^  Vg,  |Aw|  is  no  longer  a  linearly  increasing  function  of  ct.  In  fact, 
Eq.  (6)  indicates  that  Aw  will  oscillate  as  a  function  of  ct.  Consider  a  plasma  wave  of 
the  form  6kp{C,T)  =  k'ol^n/nol  cos  fcpo(C  —  Avt).  Assuming  the  laser  pulse  (with  a  pulse 
length  L  <<  Xpo  =  27r /kpo)  is  initially  centered  about  (o  =  0,  Eq.  (6)  indicates  that 
the  maximum  frequency  shift  is  given  by  (Aw,„|  =  |wpoVp«Jno/(wonoAv))  and  occurs  when 
CT  =  lApo/(2Av)l.  Furthermore,  notice  that  Awm  may  either  be  positive  or  negative, 
depending  on  the  sign  of  Av  (for  the  present  example,  Awm  >  0  for  Vp  <  Vg).  Again,  the 
assumption  that  6((,t)  is  slowly  varying  implies  |cAwTn/(2rpWo)P  <<  1. 

Physically,  the  frequency  shifts  induced  in  a  radiation  pulse  by  a  plasma  wave  may 
be  understood  as  follows.  A  plasma  wave  gives  rise  to  variations  in  the  plasma  parameter 
kp{(,T).  These  variations  lead  to  variations  in  the  “local”  phase  velocity  of  the  laser 
pulse.  Heuristically,  the  “local”  dispersion  relation  for  the  radiation  field  is  given  by 
w®«,t)  =  c^k^{(,T)  +  c^k^iCr).  For  small  ct,  this  gives  v^(C)/c^  =  1  -f  kl(,n{C)/{k^no). 
For  example,  the  local  phase  velocity  near  the  front  of  the  laser  pulse  (C  =  Co)  will  be  less 
than  the  local  phase  velocity  near  the  back  of  the  pulse  (C  =  Cl)  provided  n(Co)  <  MCl)- 
Hence,  the  individual  phase  peaks  in  the  pulse  a(z,i)  =  |a|exp(iA;;:  —  iwt)  may  move 
relative  to  one  another  (i.e.,  closer  together  for  the  present  example).  In  such  a  way  both 
the  frequency  and  the  wavenumber  of  the  radiation  pulse  will  change,  as  is  given  by  Eqs. 
(4)  and  (5). 

It  should  be  pointed  out  that  due  to  the  local  nature  of  the  frequency  shift  (the 
dependence  of  Aw  on  C))  a  laser  pulse  with  a  finite  pulse  length  L  will  develop  a  spread 


in  frequency  shifts.  That  is,  the  frequency  shift  at  the  front  of  the  pulse  may  be  different 
than  the  frequency  shift  at  the  back  of  the  pulse.  For  example,  consider  a  plasma  wave 
with  t’p  =  Vg  of  the  form  8k'^  =  fcpQ(6no/no)  sin  fcpoC  laser  pulse  centered  al^out  (  =  tt 

with  L  <  Apo/2  (see  Fig.  1).  The  difference  in  the  frequency  shift  at  the  center  of  the  pulse 
with  a  point  AC  away  is  given  by  |(Aa»(C)  —  Au>(C  +  AC))/Aa;(C)l  =  1  -  cos  AC-  This 
spread  in  frequency  sliifts  may  be  quite  significant. 

Equations  (4)  and  (5)  provide  convenient  analytic  expressions  for  the  frequency  and 
wavenumber  shifts  resulting  from  variations  in  the  plasma  parameter  k^.  However,  these 
expressions  are  only  valid  for  a  slowly  varying  envelope  which  implies  that  the 

frequency  shifts  must  be  smadl,  jAw/woP  <<  1.  To  calculate  the  behavior  of  the  vector 
potential  for  large  shifts,  |Au>/a)ol  ^  il’  is  necessary  to  retain  the  ([di^dr)  derivative 
in  the  wave  equation  for  6(C,t).  Keeping  tliis  derivative,  the  wave  equation  can  formally 
be  solved  for  the  case  Skp{C,T)  =  5fcp(C),  giving 

KCit)=  f  ‘iC'-f’(C')exp[iu;o(C' -  OKj-^oKC'.C,^)],  (9) 

Jo 

\  1/2 

where  Jo  is  the  Bessel  function  of  order  zero  and  F{C)  =  {d/d(  +  iu)o/vj)feo(C)-  in  solving 
for  6(C,t),  it  has  been  assumed  6(C  =  0)  =  0,  i.e.,  C  =  0  is  chosen  to  be  in  front  of 
the  pulse.  Analysis  of  the  above  equation  indicates  that  now  both  the  amplitude  mCji")! 
and  the  phase  evolve  as  a  function  of  r.  Determination  of  the  precise  behavior  of  b{^,T) 
requires  further  computation. 


tt(C'iCi'r)  =  I  2cV 


i: 


In  the  analysis  of  Ref.  5,  heuristic  arguments  based  on  Lorentz  transformations  of  the 
dispersion  relation  are  used  to  determine  an  expression  for  dui/dz.  The  above  analytic 
theory  may  be  used  to  derive  an  equation  for  by  operating  on  Eq.  (4)  with  d/dz  = 

d/dC-  For  the  case  of  a  plasma  wave  with  phase  velocity  Vp,  Eq.  (4)  gives 


dz 


2u;o  An 


d  d 


(10) 


which  for  Vp  =  Vg  reduces  to  the  expression  given  in  Ref.  5  when  8k^  =  k^^Snln^.  Hence, 
the  expression  given  in  Ref.  5  has  been  shown  to  be  valid  provided  i)  Vp  —  Vg,  ii)  Vp  is 
independent  of  z  and  t,  iii)  relativistic  effects  associated  with  7(2,  t)  are  neglected  and  iv) 
the  frequency  shift  remains  small,  |Au)/u;op  <<  1.  Whether  or  not  the  arguments  of  Ref. 
5  are  valid  for  large  shifts  |Au)/u)oi  >  1  remains  uncertain. 
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The  analytic  theory  presented  above  describes  how  a  laser  pulse  becomes  modified  due 
to  variations  in  the  background  plasma.  Specifically,  arbitrary  variations  in  the  plasma 
parameter  iplCiT)  leads  to  shifts  in  both  the  frequency  and  wavenumbers  indicated  by 
Eqs.  (4)  and  (5).  Provided  these  shifts  are  sufficiently  small,  the  amplitude  of  the  vector 
potential  remains  unchanged  and  propagates  forward  at  the  group  velocity  Vg.  The  above 
theory  indicates  that  changes  in  the  plasma  density  as  a  function  of  time  lead  to  shifts  in  the 
radiation  frequency.  In  particular,  the  possibility  of  upshifting  the  frequency  of  the  laser 
pulse  using  a  plasma  wave  has  been  examined  and  it  is  found  that  maximum  frequency 
shifts  result  when  Vp  =  Vg.  For  this  case,  Au>  increases  linearly  with  the  propagation 
distance  VgT  and  positive  frequency  shifts  require  phasing  the  laser  pulse  such  that  it 
is  centered  at  a  position  of  decreasing  density  or,  more  precisely,  at  a  position  where 
dSkl/dC  <  0.  The  above  theory  is  valid,  in  general,  for  small  frequency  shifts,  |Au;/a>o|^  <  < 
1.  However,  Eq.  (8)  implies  that  large  shifts  of  order  lAu>|  ~  u>o  require  propagation 
distances  of  order  cr  ~  \6no/no\~^(^^/u;pQ.  Since,  typically  (for  a  KrF  laser),  >>  1, 

this  propagation  distance  can  be  quite  large.  The  detailed  examination  of  large  shifts 
(Aw(  ~  Wo  requires  further  evaluation  of  Eq.  (9). 

This  work  was  supported  by  the  Department  of  Energy  and  the  Office  of  Naval  Re¬ 
search. 
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